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ABSTRACT 

The very young open cluster (OC) NGC 2244 in the Rosette Nebula was studied with 
field-star-decontaminated 2MASS photometry, which shows the main-sequence (MS) 
stars and an abundant pre-MS (PMS) population. Fundamental and structural pa- 
rameters were derived with colour-magnitude diagrams (CMDs), stellar radial density 
profiles (RDPs) and mass functions (MFs). Most previous studies centred NGC 2244 
close to the bright KOV star 12Monocerotis, which is not a cluster member. Instead, 
the near-IR RDP indicates a pronounced core near the 05 star HD 46150. We de- 
rive an age within 1 — 6Myr, an absorption Ay = 1.7 ± 0.2, a distance from the Sun 
d e = 1.6 ± 0.2 kpc (« 1.5 kpc outside the Solar circle), an MF slope \ = 0-91 ± 0.13 
and a total (MS+PMS) stellar mass of ~ 625 M Q . Its RDP is characterised by the 
core and cluster radii R c w 5.6' (« 2.6 pc) and -Rrdp ~ 10' (~ 4.7 pc), respectively. 
Departure from dynamical equilibrium is suggested by the abnormally large core ra- 
dius and the marked central stellar excess. We also investigate the elusive neighbour- 
ing OC NGC 2239, which is low-mass (m M s+PMS ~ 301 M Q ), young (5 ± 4Myr) 
rather absorbed (Ay = 3.4 ± 0.2), and located in the background of NGC 2244 at 
d® = 3.9 ± 0.4kpc. Its RDP follows a King-like function of R c « 0.5' k 0.5 pc and 
Rrdp « 5.0' 5.6pc. The MF slope, % = 1-24 ±0.06, is essentially Salpeter's IMF. 
NGC 2244 is probably doomed to dissolution in a few 10 7 yr. Wide-field extractions 
and field-star decontamination increase the stellar statistics and enhance both CMDs 
and RDPs, which is essential for faint and bright star clusters. 

Key words: ( Galaxy:) open clusters and associations: general; ( Galaxy:) open clus- 
ters and associations: individual: NGC 2244 and NGC 2239. 



1 INTRODUCTION 

Still in the process of emerging from the parent molecu- 
lar cloud, star clusters younger than about 5Myr usually 
present a developing main sequence (MS) and a significant 
population of pre-MS (PMS) stars. However, depending on 
the initial cluster mass, star-formation efficiency and mass 
of the more massive stars, the rapid early gas removal (from 
supernovae and massive star winds) may impart important 
changes to the original gravitational potential. One conse- 
quence of the reduced potential is that stars, especially the 
low mass ones, moving faster than the scaled-down escape 
velocity may be driven into the field. Over a time-scale of 
10 — 40 Myr, this effect can dissolve most of the very young 
star cluste rs (e.g. iGoodwin fc Bas tian 2006). Indeed, esti- 
mates (e.g. lLada fc Ladall2003l ) predict that only about 5% 



* E-mail: charles@if.ufrgs.br 
f E-mail: bica@if.ufrgs.br 



of the embedded clusters are able to dynamically evolve into 
bound open clusters (OCs). 

On observational grounds, the dramatic changes in the 
potential affecting the early cluster spatial structure should 
be reflecte d on the stellar radial dens ity profile (RDP). 
Bochum 1 (Bic a. Bonatto fe Dutrall200ct ). for instance, can 
be an example of this scenario, in which the irregular RDP 
does not follow a cluster-like profile. This suggests signifi- 
cant profile erosion or dispersion of stars from a primordial 
cluster. 

In the present paper we address the case of NGC 2244 in 
the Rosette Nebula, whi ch is also related to the Mon ocero- 
tisOB2 association (e.g. iRoman- Zufiiga fc Ladall2008D . His- 
torically, in colour-magnitude diagram (CMD) studies some 
authors centred the large-scale structure on 12 Mon, which 
is a bright foreground star of spectral type KOV. When only 
wide-field CMDs are considered, the adoption of this centre 
does not affect the results. However, as will be explored in 
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this work in the context of investigating the cluster struc- 
ture, that region is definitely at the cluster periphery. 

Based on Shanghai Observatory plates with baselines of 
34 and 87vr. lChen. de Griis fc Zhaol i|2007h derived proper 
motions (PMs) and membership probabilities for NGC 2244. 
They found mass segregation, but no velocity-mass depen- 
dence, indicating a primordial mass segregation related to 
the star-formation process. With arguments based on pub- 
lished initial mass functions (IMFs) and the measured inter- 
nal velocity dispersion of ~ 35kms -1 , they concluded that 
NGC 2244 will be dissolved on a short time-scale. 

Additionally, in the area of the Rosette Nebula, the clus- 
ter candidate NGC 2239 has been frequently included in cat- 
alogues, but hardly studied. Both NGC 2244 and NGC 2239 
are optical clusters, while the area includes numerous in- 
frare d embedded clusters in the Rosette Molecular Cloud 
(e.g. iPhelps fc Ladall 19971 ). 

This work employs 2MASS3 near-IR J, H, and K s pho- 
tometry. The 2MASS spatial and photometric uniformity al- 
low extractions of wide surrounding fields that provide high 
star-count statistics. This property makes 2MASS an excel- 
lent resource to extract photometry of a broad variety of star 
clusters, the wide field ones in particular. For this purpose 
we have been developing quantitative tools to statistically 
disentangle cluster evolutionary sequences from field stars 
in CMDs. Decontaminated CMDs have been used to inves- 
tigate the nature of cluste r candidates and derive their astro - 
physical parameters (e.g. iBica. Bonatto fc Camargol f2008). 
Basically, we apply (i) field-star decontamination to measure 
the statistical significance of the CMD morphology, which is 
fundamental to derive reddening, age, and distance from the 
Sun, and (it) colour- magnitude filters, which are essential for 
intrinsic stellar RDPs, as well as luminosity and mass func- 
tions (MFs). In particular, the use of field-star decontami- 
nation in the construction of CMDs has proved to constrain 
age and distance more than working with raw (observed) 
photo metry, especially for low-latitude OCs (|Bonatto et al.l 
12006 j ). 

2MASS can be deep for nearby young or old 
OCs. For instanc e, our group has studied the young 
OCs NGC 6611 (iBonatto. Santos Jr. fc Bical l2006h and 
NGC 4755 (|Bonatto et al.ll2006bl ). Abundant pre-MS (PMS) 
stars were seen in the w 1 Myr old NGC 6611, which is essen- 
tially embedded, and a few remaining ones in the w 14 Myr 
old NGC 4755. As ne arby older PCs we cite NGC 2 477 
jBonatto fc Bicall2005h and M67 jBonatto fc Bicall2003l ). 

In this paper we apply our set of analytical tools to the 
2MASS photometry of the stars in the area of NGC 2244 to 
derive its fundamental parameters, structure and fraction of 
MS and PMS stars. We also investigate the neighbouring 
cluster NGC 2239. 

This paper is organised as follows. In Sect. [2] we re- 
call literature data on NGC 2244. In Sect. [3] we describe 
the 2MASS photometry and compare it with the available 
optical data; we also describe the field-star decontamina- 
tion algorithm and build CMDs. In Sect. [3] we derive cluster 
fundamental parameters. In Sect. [5] we derive structural pa- 
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rameters by means of stellar RDPs. In Sect. [6] we provide 
estimates of cluster mass. In Sect. [7] we compare the struc- 
tural parameters and dynamical state of the present clusters 
with those of a sample of nearby OCs. Concluding remarks 
are given in Sect. [8] 



2 PREVIOUS WORK ON NGC 2244 

Several studies on NGC 2244, especially photometric and 
spectroscopic ones, are available in the literature. 

The WEBDAlfl database locates the cluster centre at 
a(2000) = 06 h 31 m 55 s and 5(2000) = +04°56'30", and pro- 
vides a distance from the Sun cLq = 1.45 kpc, reddening 
E(B - V) = 0.46, and an age of 7.9 Myr. 

With UBV photometry, lOgura fc Ishidal (|l98ll ) ob- 
tained E(B — V) = 0.47, a total to selective ex- 
tinction ratio Rv = 3.2, c?q = 1.42 kpc, 4 Myr of 
age, and a total mass of 5000 Mp . Wit h similar data, 
iHensberge. Pavlovski fc VerschuererJ (|200Gt ) derived an age 
of 2 Myr, R v = 3.2 ± 0.07, E(B - V) = 0.44, and d© = 
1.4 ±0.1 kpc. 

In a comprehen sive study of the Northern Monoceros 
region. Iperej |l991) found an apparent diameter of 24', an 
age within 1.45-3.63 Myr, d Q = 1.67 kpc, E(B - V) = 0.48, 
and a total mass of 770 M a for NGC 2244. 

IPark fc Sungj (|2002h found an average E(B-V) = 
0.47 ± 0.04, R v = 3.1 ± 0.2, and d Q = 1.66 kpc. By com- 
paring their photometric results with theoretical evolution 
models, they derived a main-sequence turnoff (MSTO) age 
of 1.9 Myr and a PMS age spread of about 6 Myr. The IMF 
slope calculated for the mass range 3.2 < m(Mo) < 100 is 
flat (y = -0.3 ±0.1). 

iBerghofer fc Christian! (|2002h used optical photometry 
of X-ray selected stars to estimate an isochrone age of 3 Myr, 
but significan tly younger stars are detected. 

Ill l|2005h provided updates on the nature of this young 
OC, including its central position, physical scale, and stel- 
lar population. They found substructures in NGC 2244 with 
2MASS, in particular a companion 6.6 pc west (in fact 
NGC 2239) of the NGC 2244 centre, and probably a major 
stellar aggregate resembling an arc in structure right below 
the core. Also, a disc fraction of 21 ± 3% was estimated for 
members with masses above 0.8 Mq. 

In a Chandra study of NGC 2244, IWang et all (|2008l ) 
detected over 900 X-ray sources, 77% of which having opti- 
cal or FLAMINGOS NIR counterparts. Their X-ray-selected 
population is estimated to be nearly complete bet ween 0.5 
and 3 Mq . The K-band LFs indicate a normal ISalpeterl 
l|l955l ) IMF for NGC 2244, which differs from the top-heavy 
one reported in earlier optical studies that lacked a good 
census of < 4Mq stars. The X-ray LF indicates a popula- 
tion of ~ 2000 stars with a spatial distribution strongly con- 
centrated around the central 05 star, HD 46150. The other 
early O star, HD 46223, has few companions. The cluster's 
stellar RDP shows two structures: a power-law cusp around 
HD 46150 extending for ~ 0.7 pc, surrounded by an isother- 
mal sphere reaching out to 4 pc with core radius _R C = 1.2 pc. 
This double structure, combined with the absence of mass 
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Figure 1. Top-left panel: 1° X 1° DSS B image of NGC2244 and the Rosette Nebula, centred on the coordinates of Table Q] Right: 
30' X 30' I image of the same region; circles indicate the cluster radii (Sect. O of NGC 2244 (East) and NGC 2239 (West). Several centres 
adopted for these clusters are identified in the right panel. The KOV star 12Mon is indicated by position 'B'. Other positions as indicated 
by the keys in Table[T] Bottom: K s image of NGC 2244 taken from the 2MASS Image Service focusing on the compact core within 3' X 3'. 
Same centre as in the B image. Orientation: North to the top and East to the left. 



segregation, indicates that NGC 2244 is not in dynamical 
equilibrium. The fraction of X-ray-selected members with K- 
band excesses caused by inner protoplanetary discs is 6%, 
slightly lower than the 10% disc fraction estimated from 
FLAMINGOS. The Rosette X-ray spectra of OB stars are 
soft and consistent with the standard model of small-scale 
shocks in the inner wind of a single massive star. 

Recently, I Roman- Zufiiga fc Ladal (|2008T) reviewed the 
Rosette Complex, in particular they refer to the position of 
two OCs, NGC 2244 at position 'C in the top-right panel of 
Fig-Hi an d the other at postion 'A', designated as NGC 2237. 
We call attention that their NGC 2244 actually is in the 



region of the original NGC 2239, while NGC 2237 refers to 
NGC 2244, near the position of HD 46150. 

iRoman-Zuniga fc Ladal (|2008h build a scenario where 
an expanding HII region generated by a large OB associa- 
tion interacts with a giant molecular cloud, which harbours 
a number of embedded and open clusters. 

The wealth of papers on NGC 2244 reflects the com- 
plex - and, at the same time beautiful - nature of the inter- 
play between bright massive stars, faint pre-MS stars and 
a thinning dust shroud, all embodied in a single and rela- 
tively nearby object. Fig. [T] illustrates this scenario. In the 
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Table 1. Previously adopted centres of NGC 2244 and NGC2239 
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Table Notes, (f): 12Mon as the centre of NGC 2244. Col. 7: cluster radius. 



DSSEI B image (top-left panel) NGC 2244 emerges from the 
thin dust of the Rosette central part, which is also sur- 
rounded by strong gas emission. Indeed, gas emission and 
dust absorption are nearly absent in the XDSS I image (top- 
right), and especially in the 2MASS K 3 image (bottom). 
Different centres adopted for NGC 2244 in previous (mostly 
optical) studies are indicated in the top-right panel. How- 
ever, when seen in K s , the cluster is highly concentrated 
on HD 46150 (bottom panel), suggesting a compact core. A 
sim ilar centre for NGC 2244 h a d alre ady been suggested by, 
e.g. iPerez. The fc Westerlundl ill987l ). 

The different centres are summarised in Table [T] which 
shows some confusion in the identification of the actual cen- 
tre of NGC 2244. As will be discussed in Sect. \5\ we take 
as centre the coordinates that present the maximum stellar 
density (Fig. 1 1 1 p computed within circles of 0.25' in radius, 
for MS and PMS stars taken isolately (Sect. E3J. The re- 
sulting coordinates (Table [1]) are similar to those given by 
WEBDA. The same procedure was applied to find the centre 
of NGC 2239 (Sect. [5]). 



3 NEAR-IR AND OPTICAL PHOTOMETRIES 
COMPARED 

Since the Rosette Nebula reaches about 1°, it is interesting 
to compare large-scale properties of the optical data with 
those in the near-IR. 

2MASS J, H, and K„ photometry was extracted in a 
wide circular field with VizieE0. The basic condition is that 
the extraction radius i? cxt should be large enough to al- 
low determination of the background level (Sect. [5}. We 
used i? cxt = 80' (NGC 2244) and R cxt = 30' (NGC 2239), 
which are considerably larger than the respective cluster 
radii (Sect. [5] and Tabic (2|. In the absence of significant 
differential absorption (|Bonatto fc B ica 2007a|), wide extrac- 
tion areas provide statistics for a consistent colour and mag- 
nitude characterisation of field stars. For decontamination 
purposes, comparison fields were extracted within wide rings 



located beyond the cluster radii. As photometric quality 
constraint, the 2MASS extractions were restricted to stars 
(i) brighter than the 99.9% Point Source Catalogue com- 
pleteness limii|3 in the cluster direction, and (ii) with er- 
rors in J, H, and K 3 lower than 0.1 mag. The 99.9% com- 
pleteness limits refer to field stars, and depend on Galac- 
tic coordinates. Figure [5] (panel a) shows the distribution 
of uncertainties as a function of magnitude for the stars 
in the direction of NGC 2244. The fraction of stars with 
J, H, and K s uncertainties lower than 0.05 mag is ~ 80%, 
~ 70% and ~ 60%, respectively. We employ the relations 
Aj/Av = 0.276, A H I A v = 0.176, A K JA V = 0-118 , 
and Aj = 2.76 x E(J - H) (|Dutra. Santiago fc Bicall2002l ). 
with Ry = 3.1. They stem from the extinction curve of 
ICardelli. Clayton fc Mathisl (|l989T l. 

The available B and V photometry for NGC 2244 was 
taken from SIMBAE0 within the same extraction radius as 
that used for 2MASS. As expected, the number of detected 
stars at a given radius in the optical is significantly lower 
than in the near-IR (Fig. [2] panel b). Indeed, the ratio of 
the number of stars detected in the near-IR to the optical 
NNin/Nopt increases with distance to the cluster centre, be- 
ing NMin/Nopt ~ 2 for 7? < 10' (approximately the cluster 
region) and N N i R /N op t ~ 24 for R < 80'. Except for the 
innermost region, the stellar spatial distribution detected 
with 2MASS follows a cluster-like profile (Sect. [5}, while 
the optical distribution deviates by a large amount (panel 
b). One conclusion is that analysis based on star-counts in 
a dust-rich region is more realistic in the near-IR than in 
the optical. Also, panel (b) shows that dust is thicker at 
large radii. This may have introduced biases in some of the 
previous optical studies. 

The spatial dependence of the colours towards 
NGC 2244 is examined in Fig. g] for the 2MASS (panel c) 
and optical (e) bands. The fiducial lines have been built as 
running averages of the raw (observed) data, with 10 points 
for R < 2', 100 for 2' < R < 20', and 1000 for R > 20'. 
Colours in both domains present a similar pattern, charac- 
terised by a blue core (R < 2') containing essentially the MS 



3 Extracted from the Canadian Astronomy Data Centre 
( CADC), at \http://cadcwww.dao.nrc. ca/\ 
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Figure 2. Panel (a): cumulative J, H and K s photometric error 
distribution, (b): number of stars detected by 2MASS (empty 
circles) and SIMBAD (filled squares) within a given radius; except 
for the innermost bin, the 2MASS distribution follows a cluster- 
like profile (Sect.[5]l. (c) and (d): spatial dependence of the near-IR 
colours and magnitudes, (e) and (f): same for the B and V bands. 



stars. For larger radii, foreground stars dominate the optical 
photometry, while the near-IR probes deeper regions. A sim- 
ilar effect occurs in the average magnitudes (panels d and 
f)- 

As suggested by the AR — 6' linear extractions along 
the N-S and E-W directions (Fig. [3} , a reasonable spatial 
uniformity level occurs with the 2MASS near-IR photome- 
try. Besides NGC 2244 itself, the next conspicuous bump is 
caused by NGC 2239 at » 15' to the West. These profiles 
guided the comparison field selection. 

Finally, in Fig. [3] (top panels) we show the spatial 
distribution of the stellar surface-density (a, in units of 
stars arcmin" 2 ) around NGC 2244, measured by 2MASS 
photometry. The surface density is computed in a rectan- 
gular mesh with cells of dimensions 2.5' x 2.5', with meshes 
reaching total offsets of |Aa| = |A<5| ~ 20' with respect 
to the centre (Table [TJ, in right ascension and declination. 
The respective isopleth surfaces are shown in the bottom 
panels, in which NGC 2239 shows up as a lower concen- 
tration at « 15' to the West of NGC 2244. Two cases are 
considered in Fig. [4] the observed (raw) photometry (left 
panels) and the MS + PMS stars taken separately (right) 
by means of a colour-magnitude filter (Sect. r3T3)) . Since an 
important fraction of the contaminant stars are excluded by 
the colour-magnitude filter, the surface density distribution 
of NGC 2244 (and NGC 2239) is better defined with respect 
to the surroundings. 




80 60 40 20 -20 -40 -60 -80 

AR (arcmin) 

Figure 3. Linear extractions centred on NGC 2244. The extrac- 
tions are 6' wide in both directions. The clusters are ss 15' apart 
along the EW direction. 
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Figure 4. Top panels: stellar surface-density <r(stars arcmin -2 ) 
of NGC 2244, computed for a mesh size of 2.5' X 2.5', centred on 
the coordinates in Table [T] Bottom: the corresponding isopleth 
surfaces. Left: observed (raw) photometry. Right: MS and PMS 
stars selected by means of the colour-magnitude filter (Fig. ©. 
NGC 2239 shows up at rj 15' West of NGC 2244. Act and A<5 in 
arcmin. 
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NGC2244 - Decont. NGC2239 - Decont. 
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Figure 5. Similar to Fig.[4]for the decontaminated photometries 
of NGC 2244 (left panels) and NGC 2239 (right). NGC 2239 is the 
concentration 15' to the West of the prominent NGC 2244. 

3.1 Colour-magnitude diagrams with 2MASS 
photometry 

CMDs displaying the Jx ( J — K s ) and K s x(J — K s ) colours 
built with the raw photometry of NGC 2244 are shown in 
Fig. [6] (panels a and b). The sampled region (R < 5') cor- 
responds to about half the cluster radius (Sect. [5}. When 
qualitatively compared with the CMDs extracted from the 
equal-area comparison field (panels c and d), features typ- 
ical of a very young OC emerge. A relatively vertical and 
populous MS (at 0.0 < (J — H), (J — K 3 ) < 0.3) truncated 
for stars fainter than ~ 12.5 (or mass < 2.8Mq - Sect. [4]) in 
both J and K s , stand out over the field contamination. 




1 ~ m un i is i rrrmrfKm j 

-0.3 0.1 0.5 0.9 1.3 1.7 -0.3 0.1 0.5 0.9 1.3 1.7 



(J-K s ) (J-K s ) 

Figure 6. 2MASS CMDs of NGC 2244. Top panels: observed 
photometry extracted from the R < 5' region. Middle: equal-area 
comparison field CMDs, extracted within 29.58' < R < 30', in- 
cluding the Mon Ob2 association and disc contamination. Bottom 
panels: decontaminated CMDs with the MS fitted by the 3 Myr 
Solar-metallicity Padova isochrone. PMS tracks of different ages 
are shown. The shaded polygons correspond to the MS (dark- 
gray) and PMS (light-gray) colour-magnitude filters (Sect. VS. '.it . 
The bright stars listed in SIMBAD are identified as circles; 4 of 
these are within 5 < R.(') < 10. Arrows in the bottom panels 
show the reddening vector computed for Ay = 2. 



3.2 Field decontamination 

As expected of low-Galactic latitude clusters (Table [TJ, 
the stellar surface-density in the direction of NGC 2244 
(Fig. [4} confirms that the field-star contamination, including 
Mon OB2 and disc stars, should be taken into account. Fur- 
ther confirmation is provided by the qualitative comparison 
between the CMDs extracted within the cluster and field 
(Fig. [6]l . Obviously, the field contribution should be quanti- 
fied for a better definition of the intrinsic CMD morphology. 

Although difficult, decontamination is a very impor- 
tant step in the identification and characteri sation of star 
clust ers. Most of the different approaches (e.g. iMercer et al.l 
120051 ) are based essentially on two different premises. The 
first relies on spatial variations of the star-count density, but 
does not take into account CMD evolutionary sequences. 
Alternatively, stars of an assumed cluster CMD are sub- 
tracted according to similarity of colour and magnitude with 
the stars of an equal-area comparison field CMD. Together 
with the present one, these methods are based on photo- 
metric properties only. Ideally, more robust results on mem- 
bership determination would be obtained if another inde- 



pendent parameter, such as the PM of member and com- 
parison field stars, is taken into account. However, for PM 
to be useful the cluster shou ld be relatively nearby (e.g. 
lAlessi. Moitinho fc Diasl2003l ) and/or to have been observed 
in widely-apart epochs preferentially wit h high resolution, a s 
for the globular cluster (GC) NGC 6397 (|Richer et al.ll2008l ). 
Neither condition is fully satisfied for NGC 2244, which is 
relatively distant (Sect. [4} and was observed by 2MASS in a 
single epoch. As a consequence, only about 50% of the stars 
within R = 10' of NGC 2244 have optical PM measured 
(Sect. 1^4)1. 

Our decontaminati o n algorithm is fully described in 
iBonatto fe Bical (l2007al), iBica. Bonatto fe Camargd (|2008l ) 
and Bica fe Bonattd 1 20081 ). For clarity, we provide here 
only a brief description. The algorithm measures the rel- 
ative number densities of probable field and cluster stars in 
cubic CMD cells with axes along the J magnitude and the 
(J — H) and (J — K s ) colours. It (i) divides the full range 
of CMD magnitude and colours into a 3D grid, (ii) esti- 
mates the number density of field stars in each cell based 
on the number of comparison field stars with similar mag- 
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Figure 7. Similar to Fig. [6] for the region R < 4' (panel a) of 
NGC 2239, and the equal-area comparison field (b). Panels (c)-(f): 
decontaminated CMD with different age/distance from the Sun 
solutions. Bright stars with SIMBAD optical data are indicated 
as open circles. Arrows show the reddening vector computed for 
A v = 2. 



nitude and colours as those in the cell, and (Hi) subtracts 
the expected number of field stars from each cell. Input al- 
gorithm parameters are the cell dimensions AJ = 1.0 and 
A(J — H) — A(J — K 3 ) = 0.2; the comparison fields are lo- 
cated within R = 30' - 80' (NGC 2244) and R = 20' - 30' 
(NGC 2239). The equal-area field extractions (Figs. [FJ] and 
[Jj should be considered only for qualitative comparisons. 
The decontamination itself uses the large surrounding area 
as described above. Among other statistical tests, decon- 
taminated CMDs of star clusters have integrated Ni a ^ > 1 
jBica. Bonatto fc Camargcl |200S| ; iBica fc Bonattol 120081 ). a 
condition that is fully met for the present objects, with 
JV iCT = 13.0, 6.5, respectively for NGC 2244 and NGC 2239. 



3.2.1 Decontaminated surface densities and CMDs 

We take the decontaminated surface-density distributions 
(Fig. [5j as an efficiency indicator. For the present clusters, 
the central excesses have been significantly enhanced with 
respect to the raw photometry (Fig. [4]), while the residual 
surface-density around the centre has been reduced to a min- 
imum level. By design, the decontamination depends essen- 
tially on the colour-magnitude distribution of stars located 
in different spatial regions. The fact that the decontami- 
nated surface-density presents a conspicuous excess only at 



the assumed cluster position implies significant differences 
among this region and the comparison field, both in terms 
of colour-magnitude and number of stars within the corre- 
sponding colour-magnitude bins. This meets cluster expec- 
tations, which can be characterised by a single-stellar pop- 
ulation, projected against a Galactic stellar field. 

The decontaminated CMDs are shown in the bottom 
panels of Figs. [6] and [7] As expected, essentially all contami- 
nation is removed, leaving stellar sequences typical of mildly 
reddened young OCs, with a well-developed MS and a sig- 
nificant population of PMS stars, especially in NGC 2244. 

Although in both cases the MS width is rather tight and 
appears to be dominated by photometric errors, we cannot 
exclude the possibility of differential reddening to account 
for part of the observed spread, especially towards faint 
stars. To examine this issue we show, in Figs. [6] and [7j (pan- 
els e and f), reddening vectors computed with the 2MASS 
ratios (Sect. 13. ip for a visual absorption Av = 2, approxi- 
mately the absorption derived for NGC 2244 and NGC 2239 
(Sect. [4]). Together with the decontaminated CMDs, this ex- 
periment shows that differential reddening in both clusters 
is not significant. 

We conclude that the qualitative and quantitative ex- 
pectations of the decontamination algorithm have been sat- 
isfied by the output. In both cases, the decontaminated pho- 
tometry presents a relevant excess, with respect to the sur- 
roundings, in the surface-density distribution (Fig. [SJ . In ad- 
dition, field-decontaminated CMDs extracted from the spa- 
tial regions where the excesses occur (Figs.[6]and[7]l, present 
statistically significant cluster CMDs. 

3.3 Colour-magnitude filters 

To minimise CMD noise, we apply colour-magnitude filters 
to the raw photometry to exclude stars with colours unlike 
those of the cluster sequence. The filters are wide enough to 
include cluster MS stars and the la photometric uncertain- 
tiesQ. For very young OCs such as NGC 2244 and NGC 2239, 
we also include filters to account for the PMS population. 
The colour-magnitude filters for the present OCs are shown 
in Figs. [S]and [7J 



3.4 Proper motions 

Another indication of the star cluster nature of NGC 2239 
is provided by NOMAE0 PM data taken for the stars ex- 
tracted within the same spatial region as the 2MASS data. 
However, the correspondence between NOMAD and 2MASS 
detections is not complete, with ~ 50% of the stars detected 
with 2MASS included in NOMAD, for the colour-magnitude 
filtered photometry of both NGC 2244 and NGC 2239. 

In Fig. [5] we show histograms of the right ascension 

7 Colour-magnitude filter widths should also account for forma- 
tion or dynamical evolution-related effects, such as enhanced frac- 
tions of binaries (and other multiple systems) towards the central 
parts of clusters, sin c e such systems tend to widen the MS (e.g . 
Hurley fc Toutlll998l; iKerber et al]|2002l ; iBonatto fc Bicall2007al ; 
Bonatto, Bica & Santos Jr. 2005). 

NO- 



8 http://vizier.u-stra.sbg. fr/viz-bin/VizieR?-source=I/297 . 
MAD is based on the International Celestial Reference System 
(ICRS) with origin at the solar system barycenter. 
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Figure 8. Comparative histograms of the member+field (white) 
and comparison field stars (shaded), which was scaled to match 
the projected areas. The intrinsic distributions are shown by the 
dashed lines for the extraction R < 5' of NGC 2244 (top panels) 
and R < 4' of NGC 2239 (bottom). 



(/la cos(r5)) and decimation (/xj) PM components measured 
for the member+field and field stars. We use the same clus- 
ter and comparison field extractions as those defined for the 
decontamination process (Sect. I3.2[) . The field histograms 
have been normalised to match the cluster projected area. 
Finally, the intrinsic PM distributions are obtained by sub- 
tracting the normalised field histogram from that of the 
field+members. 

Both NGC 2244 and NGC 2239 present conspicuous PM 
excesses over the field. As expected of an OC, the intrinsic 
PM distribution of NGC 2244 (Fig. [8] top panels) is essen- 
tially Gaussian. Although somewhat less defined, a similar 
conclusion applies to NGC 2239 (bottom panels). Besides, 
NGC 2244 shares essentially the same motion as the disc- 
field/MonOB2 association, which is consistent with a rel- 
atively nearby OC. NGC 2239, on the other hand, appears 
to be located at a different distance, especially because of 
the significant shift in fi a cos(<5) between member and field 
stars. 



4 CLUSTER AGE, REDDENING AND 
DISTANCE 

The field-decontaminated CMD morphologies (Sect. lB - !) can 
be used to compute cluster fundamental parameters. Both 
NGC 2244 (Fig. |6j) and NGC 2239 (Fig. present MS and 



PMS stars that can be used as constraints. We adopt solar 
metallicity isochrones because the clusters are young and 
not far from the Solar circle (see below), a r egion essentially 
occupied by [Fe/H] » 0.0 OCs (lFriellll9951 l. 

To deal with the MS we use Padova isochrones 
l|Girardi et all 120021 ) comp uted with the 2MASS J, H and 
K a filterfl The tracks of ISiess. Dufour fc Forestirul (|2000h 
are used to characterise the PMS distributions. 

We take R e = 7.2 ± 0.3 kpc (|Bica et alj|2006h as the 
Sun's distance to the Galactic centre to compute Galacto- 
centric distances, a value derived by means of the GC spatial 
distributiorQ 

Historically, different approaches have been used to ex- 
tract astrophysical parameters from isochrone fits. The sim- 
plest ones are based on a direct comparison of a set of 
isochrones with the CMD morphology, while the more so- 
phisticated ones include photometric uncertainties, bina- 
rism, and metallicit y variations. Mo st of t hese methods 
are summarised in Navlor fc Jeffries! (|2006l ). in which a 
maximum-likelihood CMD fit method is described. We cau- 
tion that, because of the 2MASS photometric uncertainties 
for the lower sequences, a more sophisticated approach for 
isochrone fitting might lead to an overinterpretation. 

For the above reasons, fits are made by eye, with the 
MS and PMS stellar distributions as constraint. We also re- 
quire that, because of the probable presence of binaries, the 
adopted (single-star) MS isochrone should be shifted some- 
what to the left of the MS fiducial line, i.e. a median line 
that takes into account the MS spread, including the photo- 
metr ic uncertainties as well (e.g. lBonatto. Bica fc Santos Jr.l 
I2005I , and references therein). In the following we discuss the 
present clusters individually. 



4.1 NGC 2244 

The decontaminated CMD morphology of NGC 2244 
(Fig. [6]l shows a nearly-vertical MS at J < 13 and 0.0 < 
(J — K s ) < 0.4, and a population of low-mass PMS stars at 
J > 13.5 and 0.9 < (J - K,) < 1.5. Taken together, these 
stellar sequences unambiguously characterise a very young 
OC. 

Allowing for photometric uncertainties, acceptable fits 
to the decontaminated MS morphology are obtained with 
any isochrone with age in the range 1 — 6Myr. The PMS 
stars in Fig. [6] are basical ly contained within the 0.2 My r 
and 6 Myr PMS isochrones jsiess. Dufour fc ForestinfeOOOf ) . 
thus implying a similar age range as the MS. Accordingly, 
we take the 3 Myr isochrone as representative solution. 

With the adopted solution, the fundamental parameters 
of NGC 2244 are a near-IR reddening E(J - H) = 0.17 ± 
0.02 (or E(B - V) = 0.54±0.06 and A v = 1.7±0.2), an ob- 
served and absolute distance moduli (m — M)j = 11.5±0.2 
and (m — M)o = 11.03 ± 0.21, respectively, and a distance 



9 http://stev.oapd.inaf.it/cgi-bin/cmd These isochrones 
arc very similar to the Johnson-Kron-Cousins ones (e.g. 
iBessel fc Brettj|l988h. with diff erences of at most 0.01 in colour 
llBonatto. Bica fc Girardill2004) . 

10 Other recent studies gave s i milar results, e.g. Rq = 
7.2 ± 0. 9 kpc [lEisenhauer etaTI 120031) , Rq = 7.62 ± 
0.32 kpc llEisenhauer et al.l 12005^ and _R = 7.52 ± 0.10 kpc 
dNishivama et al]|2006h . with different approaches. 
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Figure 9. Colour-colour diagrams with the decontaminated pho- 
tometry of NGC2244 (top-left panel) and NGC 2239 (bottom - 
left). The 3Myr and 5Myr ISiess, Dufour fc Forestiiiil (|2000|) 
isochrones,set with the derived reddening values, are used to show 
the PMS and MS sequences simultaneously. Reddening vectors for 
A v = 1.7, 3.4 are shown for NGC 2244 and NGC 2239, respec- 
tively. Right panels: evolutionary tracks of PMS stars of differ- 
ent masses are superimposed on the decontaminated CMDs of 
NGC 2244 (top) and NGC 2239 (bottom). 



from the Sun d© = 1.6±0.2 kpc. Thus, for R e = 7.2 kpc, the 
Galactocentric distance of NGC 2244 is R G c = 8.7T0.2 kpc, 
which puts it ps 1.5 kpc outside the Solar circle. This solu- 
tion is shown in Fig. [6] (panels e and f). The age spread 
indicates a non-instantaneous star formatio n process, sim- 
ilar to what was foun d for, e.g. NGC 4755 llBonatto et all 
l2006bl ) and NGC 6611 l|Bonatto. Santos Jr. fc Bicall2006l ). 



4.2 NGC 2239 

Although somewhat less-populated than NGC 2244, the de- 
contaminated CMD of NGC 2239 also presents MS and PMS 
sequences (Fig. [7]), which suggests a similar age as that of 
NGC 2244. However, this less-constrained CMD admits al- 
ternative age/distance solutions. 

Thus, before the CMD fit, we compute the distance 
of the central bright stars projected on NGC 2239. We 
found 2 bright stars in SIMBAD in common with the 
2MASS detections. These are the A 2 star GSC 00154-01659 
(B = 13.32, V = 12.90, J = 11.336) and the B5 star 
GSC 00154-02384 (B = 13.42, V = 13.21, J = 12.014). 
With the typical col o urs a nd magnitudes of stars given in 
iBinnev, fc Merrifieldl i|l998h . the A 2 star would be located 
at d© = 1.4 kpc, while the B 5 would be at d© = 4.3 kpc. 



Together with young MS and PMS isochrones, both dis- 
tances are used as starting point to search for fundamen- 
tal parameters. The possible solutions are summarised in 
Fig. [7] for both the far (left panels) and near (right) dis- 
tances. We consider the MS ages of 1 Myr (middle panels) 
and 5 Myr (bottom). The best fit for these ages was for 
d© m 3.9 ± 0.4kpc and d© ~ 1.7 ± 0.3 kpc. Within uncer- 
tainties, they are consistent with those derived for the bright 
stars. The MS+PMS solutions for d© « 3.9 kpc (Fig. [7] pan- 
els c and f) appear to account for essentially all of the de- 
contaminated CMD sequences. Irrespective of age, the near 
solutions, on the other hand, fail to explain several faint 
stars below the MS. 

Considering all available morphological and photomet- 
ric properties, the best overall CMD fit (Fig. [7] panel e) 
corresponds to an age 5 ± 4 Myr, E(J - H) = 0.34 ± 0.02, 
E{B — V) = 1.10 ± 0.06, Av = 3.4 ± 0.2, {m-M)j = 
13.9 ± 0.2, (m - M)o = 12.95 ± 0.21, d Q = 3.9 ± 0.4kpc, 
and J?gc = 10.8 ± 0.3 kpc, thus » 3.6 kpc outside the So- 
lar circle. Thus, the B5 star GSC 00154-01659 is probably 
a cluster member, while the A 2 star GSC 00154-01659 is 
probably in the foreground. We conclude that NGC 2239 
lies « 2.3 kpc in the background of NGC 2244. 



4.3 Colour-colour diagrams 

When transposed to the near-IR colour-colour diagrams 
(J — K s ) x (H — K s ), the age and reddening solutions of 
NGC 2244 and NGC 2239 derived above consistently match 
the field-star decontaminated photometry of these OCs 
(Fig. |9j le ft panels). Since they inc l ude PM S stars, we use 
tracks of ISiess. Dufour fc Forestinil (2000) to characterise 
the age of NGC 2244 (~ 3 Myr) and NGC 2239 (~ 5 Myr). 
MS stars lie on the blue side of the diagram. Interestingly, 
only a low fraction of the stars of these OCs appears to 
be very reddened. In both cases, the distribution suggests a 
low fraction of d isc stars, consisten t with the estimates in 
the range 6-10% jWang et alj|2008h . 



4.4 The fraction of PMS and MS stars 

Right after formation, most of a cluster's mass is stored 
in the PMS stars that, eventually, shed the dust lay- 
ers and emerge into the MS (Sect. [1}. Thus, the num- 
ber of MS and PMS stars evolve in opposite directions 
with clu ster age, till all star s are in the MS after about 
30 Myr l|Bonatto et al.l l2006bl . and references therein) . In- 
deed, NGC 2244 and NGC 2239 present different fractions 
of MS and PMS stars (Figs. [6]and[7)l, which is consistent 
with the different ages. 

To further explore this issue we compute the ratio of 
the number of PMS to MS stars fpMS/MS = npMs/nMS- 
Then we examine the age dependence of fpus/US f° r 
the very young OCs, located at a similar distance as 
NGC 2244, studied by our group with the same methods 
as those employed in NGC 2244 and NGC 2239. These con- 
ditions are satisfied by NGC661]E1 (1-3 ± 0.3 Myr, d© w 



For consistency with the remaining clusters, we counted the 
number of PMS stars brighter tha n J = 16 in NGC 6611, since in 
iBonatto. Santos Jr. fc Bical j2006l) we restricted them to J ^ 15. 
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Figure 10. The ratio of the number of PMS to MS stars fol- 
lows an exponential-decay function (dashed line) with cluster age: 
n PMs/ n MS exp(— age/r), with r = 5.6 ± 1.6 Myr. Fit uncer- 
tainties (lex) are within the shaded region. 



1.8 kpc— IBonatto, Santos Jr. fc Bicall2006h. Bochum 1 (9± 



Bica. Bonatto fc Dutrall2008l). and 



1.8 kpc — IBonatto et al.l 



3 Myr, d Q « 2.0 kpc 
NGC4755 (14 ± 4 Myr, d© 
l2006bh . 

Despite the considerable uncertainties in age and the 
relatively small number of clusters, the fpMS/MS ratios 
shown in Fig. [TD] appear to be consistent with the expected 
trend with cluster age. Indeed, we found that the ratios fol- 
low the exponential-decay function fpMS/MS K e~^ age ^ T \ 
with the time-scale r = 5.6 ± 1.6 Myr. Thus, it would take 
less than about 30 Myr for a cluster to retain only m 5% 
(a typical Poisson fluctuation) of the original PMS stars. 
It would be interesting to check this trend - and the above 
time-scale - for a statistically significant sample of young 
clusters, so that other parameters, such as cluster mass, and 
environmental effects can be considered as well. 



5 CLUSTER STRUCTURE 

We use the projected stellar RDPs, defined as the stellar 
number density around the cluster centre, to derive struc- 
tural parameters. To minimise noise, we work with colour- 
magnitude filtered photometry to isolate the MS and PMS 
stars, which enhances the RDP contrast rel ative to the back- 
groun d, especially in crowded fields (e.g. IBonatto fc Bical 
12007 j ). However, field stars with colours similar to those 
of the cluster are expected to remain inside the colour- 



magnitude filter, affecting the intrinsic RDP in a way that 
depends on the relative densities of field and cluster stars. 
The contribution of the residual contamination to the ob- 
served RDP is statistically evaluated by means of its exten- 
sion into the field. 

Rings of increasing width with distance from the cluster 
centre are built to avoid oversampling near the centre and 
undersampling at large radii. The set of ring widths used is 
AR = 0.25, 0.5, 1.0, 2.0, and 5', respectively for 0' sC R < 
0.5', 0.5' < R < 2', 2' < R < 5', 5' < R < 20', and R ^ 
20'. The residual background level of each RDP corresponds 
to the average number-density of filtered field stars. The R 
coordinate (and uncertainty) of each ring corresponds to the 
average position and standard deviation of the stars inside 
the ring. 

The colour-magnitude filtered RDPs of the clusters are 
shown in Fig. QT] As expected, minimisation of the number 
of non-cluster stars by the colour-magnitude filter resulted 
in RDPs with a high contrast relative to the background. 
For NGC2244 we also show the RDPs built with the MS 
and PMS stars separately (left panels). Interestingly, while 
the MS RDP (panel b) has a conspicuous density excess 
for R rs 0.15', the PMS stars (panel c) are found only for 
R > 0.4'. The presence of NGC 2239 causes a bump in the 
MS RDP of NGC 2244 (panel b). Similarly, NGC 2244 shows 
up in the RDP of NGC 2239 (d). Fig.[TT]also shows the RDP 
produced with the star 12 Mon as centre (panel d). It is clear 
that 12 Mon cannot be the centre of NGC 2244. 

Most star clusters have RDPs that follow a well-defined 
analytical profile e.g. , the empir ical, single mass, m odified 
isothermal spheres of lKind (1 19661) and Wilson! lU975|fp|. an d 
the power law with a core of lElson. Fall fc Freemarl l|l987h . 
Each function is characterised by a different set of parame- 
ters that are related to cluster structure. For simplicity and 
considering the error bars of the RDPs in Fig. 111! we adopt 
the function a(R) = ab g + co/(l + (R/R c ) 2 ), where o\,g is 
the residual background density, gq is the central density of 
stars, and 7? c is the core r adius. It is similar to the function 
introduced bv lKind (|l962h to describe the surface brightness 
profiles in the central parts of GCs. To minimise degrees of 
freedom, oo and R c are obtained from the fit, while o^g is 
measured in the field. The RDP bins corresponding to the 
neighbouring clusters were ignored in the fit. The best-fit 
solutions are shown in Fig. 1111 and the parameters are given 
in Table [2] For absolute comparison with other clusters, Ta- 
ble [2] gives parameters in absolute units. 

Within uncertainties, the adopted King-like func- 
tion describes well the colour-magnitude filtered RDP of 
NGC 2239 (panel d) along the full radius range. The same 
applies only for R > T for the RDP of NGC 2244. The 
innermost bin in the MS (and to a lesser degree to the 
MS+PMS) RDP (panel b) presents a several a excess over 
the fit. This RDP cusp basically corresponds to the de- 
tached grouping of stars (with a diameter of « 0.25') around 
HD 46150, seen in Fig. [1] (bottom-right panel). Ou r inner 
RDP shape agrees with that derived by I Wang et al.l (|2008l ) 
with FLAMINGOS. In old star clusters, such a central 
RDP excess can be attributed to a post-core collapse, like 



12 It assumes a pre-defined stellar di stribution fu nction and pro- 
duces more extended envelopes than King (1966). 



The very young OCs NGC22U and NGC2239 11 



Table 2. Derived cluster structural parameters 
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(pc) 
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(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


NGC 2244t 


1.53 ±0.02 


3.87 ±0.65 


5.6 ±0.8 


10.0 ± 2.0 


3.5 ±0.6 


0.466 


7.0 ±0.4 


17.8 ±3.0 


2.6 ±0.4 


4.7 ±0.9 


NGC 2244* 


0.12 ±0.01 


0.95 ±0.43 


1.2 ±0.9 


8.0 ± 1.0 


9.2±4.1 


0.466 


0.5 ±0.1 


4.4 ± 1.9 


0.6 ±0.4 


3.7 ±0.5 


NGC 2244* 


1.41 ±0.09 


3.61 ±0.61 


5.7 ±0.8 


12.0 ± 2.0 


3.6 ±0.6 


0.466 


6.5 ±0.4 


16.6 ±2.8 


2.7 ± 0.4 


5.6 ±0.9 


NGC 2239 1 " 


1.90 ±0.05 


12.77 ±4.77 


0.5 ±0.1 


5.0 ± 1.0 


7.7 ±2.9 


1.127 


1.5 ±0.1 


10.0 ±3.7 


0.5 ±0.1 


5.6 ± 1.1 



Table Notes. RDPs built considering separately the MS±PMS (f), MS (J) and PMS (*). Core (i? c ) and cluster (-Rrdp) radii are given in 
angular and absolute units. Col. 6: cluster/background density contrast parameter (<5 C = 1 ±o"o/ofc s ), measured in the colour-magnitude 
filtered RDPs. Col. 7: arcmin to parsec scale. 



10 




10 



0.1 



10 



NGC2244 - MS + PMS 

(a) 



NGC2244 - MS 
(b) 




NGC2239 - MS + PMS 
(d) 



12 MON 



1 10 

R (arcmin) 



0.1 



1 10 
R (arcmin) 



Figure 11. Stellar RDPs built with colour-magnitude filtered 
photometry. Solid line: best-fit King-like profile. Horizontal 
shaded polygon: background. Shaded regions: Icr King fit uncer- 
tainty. Note the central density excess in the RDP of NGC 2244 
in panels (a) and (b). The RDP with 12 Mon as centre in shown 
in (d). 



those detected in some GCs fe.g. lTrager. King fc Diorgovskil 
1995). It has been det ected as well in Gyr -class OCs, 
such as, e.g. NGC 3960 (|Bonatto fc Bical 120061 ) and LK 10 
|Bonatto fe Bical 120091 ). Another very young cluster har- 
bouring such a de tached core producing an RD P central 
cusp is NGC 6823 l|Bica. Bonatto fc Dutral 120081 ). Clusters 
are not expected to dynamically evolve into a post-core col- 
lapse on short time-scales, and the cusp must have been 
caused by star-forming effects. Th e compact core within 
the eroded profile of Bochum 1 (|Bica. Bonatto fc Dutral 



2008) can be a long-lived structure in young clusters. Con- 
sequently, this central cusp in such a young cluster as 
NGC 2244 suggests a significant deviation from dynamical 
equilibrium (Sect. [7]). 

We also estimate the cluster radius (-Rrdp) by visually 
comparing the cluster RDP and background levels, i.e. -Rrdp 
is the distance from the cluster centre where both are sta- 
tistically indistinguishable (e.g. iBonatto fc B ica 2005|, and 
references therein). Most of the cluster stars are contained 
within -Rrdp, which should not be mistaken for the tidal 
radiuf^l. The cluster radii of the present objects are given 
in angular and absolute scales (Table [2}. 

The density contrast parameter S c — 1 + crn / o~b g , which 
is relatively high (3.5 < S c < 9.2) for the present RDPs, 
is also given in Table [2] Since 8 C is measured in colour- 
magnitude-filtered (lower noise) RDPs, it is usually higher 
than the visual contrast produced by images (e.g. Fig. [T}. 

Taken at face value, the core radius of NGC 2244 (for 
the MS±PMS stars) R c ~ 2.6 pc would put it beyond the 
median value of the di stribution deriv e d for a sample of rela- 
tively nearby OCs bv lPiskunov et al.l l|2007l ). NGC 2239, on 
the other hand, falls on the low--R c tail. Besides, assumin g 
the relation tidal radius ~ 2 x -Rrdp l|Bonatto fc B ica 2005), 
both clusters fall around the median value of the tidal radius 
distribution. 



6 CLUSTER MASS 

Both clusters clearly present distinct populations of MS and 
PMS stars (Figs. |S]and[7|. As the first step to estimate 
the cluster masses we build the luminosity functions (LFs) 
in the K s band for the MS and PMS stars separately, by 
means of the respective colour-magnitude filters (Sect. |3"3)) . 
We show them in Fig. 1121 where the similar age, different 
distances and number of members are reflected, especially 
on the different MS and PMS cutoffs. In both cases the 
PMS LFs present the expected steep increase towards faint 



13 Tidal radii are derived from, e.g. t he 3-parameter King-profile 
fit to RDPs llBonatto fc Bical l200c^ . which requires large sur- 
rounding fields and adequate errors. For instance, in populous 
and relatively high Galactic latitude OCs such as M 67, NGC 188, 
and NGC 2477, the tid al radii are a factor ~ 2 larger than -Rrdp 
dBonatto fc Bicall2005h . 
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magnitudes (low-mass stars) , which confirms that PMS stars 
are an important fraction of the members. 

For a more objective investigation on the stellar mass 
distribution we build the MFs (</>(m) = ^) for the cur- 
rent MS stars that, in turn, can be used to compute the 
mass stored in stars. Similarly to the RDPs (Sect. [SJl, we 
work with colour-magnitude filtered photometry to minimise 
noise. First we build the LF independently for each 2MASS 
band, both for the cluster region (R < Rrdp) and compari- 
son field. The intrinsic LFs are obtained by subtracting the 
respective (equal-area) comparison field LF from that of the 
cluster. The intrinsic LFs are transformed into MFs with the 
mass-luminosity relations obtained from the corresponding 
age and distance from the Sun solutions (Sect.|4)|. The final 
MF is produced by combining the J, H and K s MFs into 
a single MF. Further details on MF construction are given 
in lBica. Bonatto fe Blumberd (120061 ). The effective MS stel- 
lar mass ranges are (4.6 ± 2.2) s; m(M Q ) < 60 (NGC2244) 
and (3.2 ± 0.5) < m(M©) < 14 (NGC2239). As Fig. [jj 
(bottom panels) shows, the MS MFs are rather smooth and 
present different lower and upper masses, which reflects the 
lower distance and younger age of NGC 2244 with respect 
to NGC 2239. 

Since PMS stars are abundant in both clusters, it is im- 
portant to build their MF as wel l. In Fig. [9] (right panels) 
we sh ow the evolutionary tracks ( Siess. Dufour fc Foresting 
2000) of PMS stars of different masses superimposed on the 
decontaminated CMDs of NGC 2244 and NGC 2239. It is 
clear, especially for NGC 2244, that PMS stars less massive 
than IMq are the most abundant component. Similarly to 
the MS, the PMS MFs are built with the number of PMS 
stars among any two tracks in the cluster region and compar- 
ison field. Finally, we add the MS and PMS MFs to produce 
the total MF of each cluster (Fig. [12]). 

The number of MS (ums) an d PMS (npA/s) mem- 
bers in NGC 2244 (for R ^ -Rrdp) are derived by counting 
the stars in the background-subtracted colour-magnitude fil- 
tered photometry. We apply the same approach as above 
to compute the PMS mass. There are ums = 26 ± 3 and 
npias = 301 ± 60 stars; the corresponding mass values are 
niMS = 389 ± 44 M and m PM s = 236 ± 46 M (com- 
puted assuming the average mass between any two evolu- 
tionary tracks in Fig. [5J. Thus, the total stellar mass of 
NGC 2244 is mus+PMS ~ 625 M , wh ich agrees with the 
770 Mq mass estimated bv lPered (|l99ll ). We note that this 
value is about 10% of the mass estimated bv lOgura fc Ishidal 
(1981) for NGC 2244. However, this difference may arise 
from the present detailed analysis - especially the decon- 
tamination and the separation of MS and PMS stars in 
the construction of the cluster MF. The same analysis ap- 
plied to NGC 2239 yields n M .s = 26 ± 3, n PM .s = 70 ± 11, 
m MS = 141 ± 16 Mq, niPMs = 160 ± 32 Mq , and the total 
stellar mass ttims+pms ~ 301 Mq, about half the mass of 
NGC 2244. 

Considering the MS stars isolately, the MFs can be 



well represented by the function 



-(i+x) 



with the 



slopes x = 0.24 ± 0.09 and x = 0.76 ± 0.08, respectively for 
NGC 2244 and NGC 2239. Both values are flatter than the 
X = 1.35 of lSalpeterl (|l955l ^ initial mass funct ion (IMF). A 
flat M F slope was also found for NGC 2244 bv lPark fc Sund 
(2002). However, when the MS and PMS stars are taken 
together, the slopes become steeper, x = 0.91 ± 0.13 and 
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Figure 12. Top panels: i^s-luminosity functions of the MS (filled 
circles) and PMS (empty circles) stars. Bottom: mass functions 
(for the J, H and K s bands combined) of the MS (filled cir- 
cles) and MS+PMS (empty squares). Fits of the function <f>(m) oc 
m -(l+x) ar e shown by the dashed (MS) and dotted (MS+PMS) 
lines. 



X = 1.24 ± 0.06. While within u ncertainties the total MF 
of NGC 2239 is comparable to the lSalpeterl (|l955l ) IMF, the 
MF of NGC 2244 remains so mewhat flatter, bu t still consis- 
tent with the conclusions of IWang" et all l|2008h . 



7 DISCUSSION 

Constrained by isochrone fits (Sect. we could derive 
fundamental and structural parameters of the young OCs 
NGC 2244 and NGC 2239, part of them for the first time. 
We use them to compare some of their properties with those 
of well-studied OCs. 



7.1 Diagnostic diagrams 

We further investigate the nature of NGC 2244 and 
NGC 2239 with diagrams that examine relations among 
astrophysical parameters of OC s in different e nviron - 
ments. They were introduced by iBonatto fc Bical ([2005) . 
As reference sample w e use some br i ght ne arby OCs 
(jBonatto fc Bical 120051 : IBonatto etail l2006bl ). and a 
group of O Cs projected towards t he central parts of 
the Galaxy (|Bonatto fc Bical l2007al ). Also included are 
the young OCs NGC 6611 wit h the age ~ 1.3 Myr 
l|Bonatto. Santos Jr. fc Bicall2006l). NGC 6823 with ~ 4 Myr 



and Bochum 1 with ~ 9 Myr l|Bica. Bonatto fc Dutrall200; 
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Figure 13. Diagrams dealing with astrophysical parameters of 
OCs. Gray-shaded circles: reference OCs. The young star clusters 
NGC 6611, NGC 6823 and Bochum 1 are indicated for comparison 
purposes. Analytical relations in panels (c) and (e) are disc ussed 
in the text. Dotted line in panel (f) shows [Salpeterl jl955Tl IMF 
slope x = 1.35. 



NGC 6611 and NGC 6823 serve as comparison with gravita- 
tionally bound objects of similar age, while Bochum 1 is a 
star cluster fossil remain that might be dynamically evolving 
into an OB association. The full sample of comparison OCs 
is characterised by ages in the range ~ 1.3 Myr to ~ 7Gyr, 
and Galactocentric distances within 5.8 < RGc(kpc) < 8.1. 
Their parameters have been obtained following the same 
prescriptions as those for NGC 2244 and NGC 2239. 

The diagrams are shown in Fig. 1131 where panels (a) and 
(b) examine the dependence of cluster (Rrdp) and core (_R C ) 
radii on cluster age, respectively. Most of the small-radius 
OCs (especially in -Rrdp) occur at an age ~ 0.5 — 1 Gyr, the 
typical tim e-scale of OC disruption processes near the Solar 
circle (e.g. iBergond. Leon fc Guil bcrt 2001; lLamers et al.l 
|2005| ). Both NGC 2244 and NGC 2239 present a cluster ra- 
dius comparable to that of NGC 6611 (and in general equiv- 
alent to other young OCs). The same applies to the core 
radius of NGC 2239. NGC 2244, on the other hand, has an 
R c too large when compared to the reference OCs. 

Core and cluster radii of the reference OCs follow the 
relation R RD p = (8.9 ± 0.3) x i?^ 10±01) (panel cQ sug- 
gesting a similar scaling for both kinds of radii. While 
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iNilakshi. Pandev fc Mohanl d2002^ . ISharma etUI (|2006h . 
and lMacieiewski fc Niedzielskil J2007T) ■ 



NGC 2239 fits tightly in that relation, NGC 2244 deviates 
again probably because of the exceeding cluster radius. A 
dependence of OC size on Galactocent ric dist a nce is sug- 
gested by pa n el (d) , as discussed by iLvngal (Il982l ) and 
iTadross et all l|2002l ). While NGC 2244 follows the trend, 
NGC 2239 deviates somewhat. This relation may be partly 
primordial, in the sense that the high molecular gas density 
in ce ntral Galactic regions may have produced small clusters 
(e.g. Ivan den Bergh. Morbev fc Pazderlll99lh . After forma- 
tion, mass loss due to stellar and dynamical evolution (e.g. 
mass segregation and evaporation), together with tidal in- 
teractions with the Galactic potential and giant molecular 
clouds, also contribute to the depletion of star clusters, es- 
pecially the low-mass and centrally located ones. 

When the mass-density r adial distributi o n fol- 
lows a King-like pr o file (e.g. iBonatto fc Bical l2007bl : 
IBonatto fc Bical 120081 ; IBonatto. Bica fc Santos Jr.l |2008h . 
the cluster mass inside -Rrdp can be computed as a function 
of the core radius (R c ) and the central mass-surface density 
(umo), M c i us = 7T R^o-Moln [l + (Rrdp/Rc) 2 ] • With the 
above relation (panel c) between R c and -Rrdp, this equa- 
tion becomes M c i us ~ 13.8ctmoRc- The observed relation 
of core radius and cluster mass is examined in panel (e). 
The reference OCs, together with NGC 2239 are contained 
within King-like distributions with central mass densities 
within 30 < ctmo (Mq pc" 2 ) < 600. NGC 2244 appears to 
have too big a core radius for its total mass. 

Finally, when the total (MS+PMS) MF slope is con- 
sidered (panel f), NGC 2239 and especially NGC 2244, ap- 
pear to have MFs flatter than those of similarly young OCs. 
On the other hand, their slopes are equivalent to those de- 
rived for some old OCs in the reference sample. In gen- 
eral, fiat MFs reflect advanced dynamical evolution (e.g. 
IBonatto fc Bicall2005l ). 

What follows from the above analysis is that while 
NGC 2239 appears to be characterised by parameters of a 
typical young OC, NGC 2244, on the other hand, exhibits 
evidence of a system far from dynamical e quilib rium, which 
agrees with IChen. de Griis fc Zhaol (|2007l ) and lWang et al.l 
i|2008h . 



8 SUMMARY AND CONCLUSIONS 

In the present paper we employ the wide-field and near-IR 
depth provided by 2MASS to focus on the Rosette Nebula 
cluster NGC 2244 and the nearby projected OC NGC 2239. 
Our approach relies essentially on field-star decontaminated 
2MASS photometry, which enhances cluster CMD evolu- 
tionary sequences and stellar radial density profiles, produc- 
ing more constrained fundamental and structural parame- 
ters. 

Previous studies were mostly based on optical photom- 
etry and/or near-IR with small angular fields. However, 
2MASS can still provide additional insight (Sect. [3}- The 
set of tools developed by our group allowed to unambigu- 
ously isolate MS and PMS stars that, in turn, resulted in 
well-defined CMDs, RDPs and mass functions. In addition, 
we explore proper motion properties to investigate the other 
cluster in the area, NGC 2239. 

Taken together, the (decontaminated) MS and PMS se- 
quences of NGC 2244 provided an age range 1 — 6 Myr, an 
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absorption Av = 1.7 ± 0.2, and a distance from the Sun 
do = 1.6 ± 0.2 kpc (~ 1.5 kpc outside the Solar circle). 
These parameters are consistent with most of the previ- 
ous estimates. The (MS+PMS ) MF slope y = 0.91 ± 0.13 
is somewhat flatter than the ISalpeterl (|l955l ) IMF. With 
a total (MS+PMS) stellar mass of itims+pms ~ 625 M 
derived in the present work, NGC 2244 is not as massive 
as previously estimated (Sect. [2}. The King- like profile fit 
to the (MS+PMS) stellar RDP was obtained with a core 
radius R c pb 5.6' pb 2.6 pc; the corresponding cluster ra- 
dius is -Rrdp ~ 10' ~ 4.7 pc. Compared to a set of well- 
known OCs, the core radius of NGC 2244 appears to be 
abnormally big, which puts it at unusual loci in the struc- 
tural/dynamical diagnostic-diagrams (Sect.[7T}. NGC 2244 
has a central cusp that cannot be fitted by e.g., a King's 
law. This cusp is probably due to a star-formation effect, 
and not the product of dynamical evolution. We conclude 
that NGC 2244 is not in dynamica l eq uilibrium, consisten t 
with lChen. de Griis fe Zhaol (|2007h and I Wang etail (|200Sh . 

NGC 2239 is a low-mass {iums+pms ~ 301 M Q ), young 
(5 ± 4Myr) and somewhat more absorbed (Av = 3.4 ± 0.2) 
OC, at d© = 3.9 ± 0.4 kpc, thus in the background of 
NGC 2244. Structurally, its RDP can be represented by a 
King-like profile with R c ps 0.5' ~ 0.5 pc and -Rrdp ~ 5.0' ~ 
5.6 pc. With % = 1-24 ± 0.06, its composite MS+PMS MF 
slope is essentially Salpeter's IMF. These parameters char- 
acterise an average young OC, as compared to the reference 
nearby OCs fSect. l7TTj) . 

While NGC 2239 is a normal young OC with MS 
and PMS stars distributed according to a cluster RDP, 
NGC 2244 appears to be ano ther example, like Bochum 1 
|Bica. Bonatto fe Dutral [2008) , of an open cluster doomed 
to dissolution in a few 10 7 yr. The present work shows the 
importance of field-star decontamination and wide-field ex- 
tractions to get the best stellar statistics and to produce 
high-quality CMDs and RDPs. 
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